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ABSTRACT 

The  thesis  develops  a  theoretical  model  for  predicting 
ship  interactions  that  should  be  useful  in  a  ship  trajec- 
tory analysis.   The  method  assumes  that  when  ships  are 
close  together,  the  potential  field  is  the  main  source  of 
the  interaction  forces.   The  flows  around  the  ships  are 
represented  by  using  half  bodies  of  revolution  generated 
with  axial  singularity  distributions.   A  second  axial 
distribution  of  doublets  were  included  to  correct  for 
cross  flows  and  thus  maintain  the  rigid  body  boundary 
condition.   Following  the  determination  of  the  fluid  flows, 
the  forces  were  calculated  using  Lagally's  theorem  as  ex- 
panded by  W.  E.  Cummins. 

The  principal  result  was  a  computerized  model  to  pre- 
dict forces  resulting  from  steady,  parallel  motion  of 
bodies  of  revolution  moving  in  an  infinite  ideal  fluid. 
The  program  is  capable  of  handling  ship  geometries  in  rec- 
tilinear motion.   It  was  found  that  the  theoretical  model 
produced  similar  results  when  compared  with  model  tests 
conducted  by  D,  W.  Taylor,  hov/ever,  those  tests  were  re- 
ported as  being  crude.   This  emphasizes  the  need  to  obtain 
a  set  of  more  valid  empirical  data. 
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CHAPTER  I 
INTRODUCTION 

1. 1   General 

Interactions  between  ships  operating  in  proximity  can 
be  thought  of  as  a  modern  day  phenomena.   There  are  no  knov/n 
references  to  it  in  the  literature  until  after  the  emer- 
gence of  steam  driven  ships.   Along  with  the  increase  in 
shipping  in  shallow  and  restricted  waters,  this  innovation 
contributed  to  the  higher  frequency  of  collisions.   These 
collisions  provided  the  operator  with  his  initial  encounters 
with  ship  interaction  forces.   It  was  many  years,  however, 
before  the  courts  recognized  the  existence  of  this  strange 
force  which  contributed  to  many  of  the  collisions.   The 
acceptance  of  these  forces  came  only  after  many  eyewitness 
accounts  described  strange  ship  behaviors  and  it  was  con- 
cluded that  there  must  have  been  a  force  that  was  not 
present  during  normal  operations. 

With  the  increased  density  of  shipping  in  restricted 
waters,  especially  for  channels,  it  is  not  surprising  that 
collisions  involving  interactions  were  first  noticed  for 
the  shallow  water  occurrences.   The  effect  of  shallow  water 
has  persisted  as  being  important  in  describing  the  inter- 
action phenomena.   Shallow  water  does  influence  the  flow 
around  the  ship,  however,  for  depths  greater  than  twice  the 

draft  of  the  ship,  the  effect  of  the  ocean  bottom  is 
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expected  to  be  small.   This  is  based  upon  the  interaction 
between  ships  in  deep  water,  forces  have  been  experimentally 
determined  as  the  ship  separation  was  increased, 
^^     With  the  increased  speeds  and  maneuverability  intro- 
duced by  steam  propulsion,  interactions  other  than  with 
ships  were  first  noticed  and  recorded.   These  involved  the 
forces  experienced  by  the  position  of  the  ship  in  the 
channel.   From  these  forces,  rules  of  thumb  were  formed  such 
as  the  ship  "smells  the  bank".   This  particular  rule  indi- 
cates the  mariner  realized  that  the  ship  interacts  with  the 
channel  bank  so  as  to  sheer  away  from  it.   It  is  in  this 
area  of  ship  maneuverability  in  restricted  water  that  a 
great  many  empirical  studies  have  been  made.   With  the  re- 
sults of  Bernoulli,  the  mechanism  could  be  explained.   Even 
today  the  operators  aware  of  the  Venturi  effect  tend  to 
associated  all  interactions  and  explain  them  consistent 
with  this  law.   This  may  give  a  rough  guide,  but  does  not 
contribute  to  the  understanding  of  the  process.   It  may  at 
times  even  lead  an  investigator  to  the  v/rong  conclusions. 
The  forces  evolving  from  ship  interactions  has  long 
been  a  controversial  subject.   The  main  interest  has  stemmed 
from  the  cases  which  result  from  collisions.   The  literature 
is  sprinkled  with  accounts  of  major  collisions.   Included 
are  discussions  where  important  figures  of  the  day  give 
conflicting  explanations  of  the  ships'  behavior.   The  first 


scientific  investigations  involving  model  tests  occurred 
in  1909  when  D.  W.  Taylor  presented  his  findings  to  the 
Society  of  Naval  Architects  and  Marine  Engineers.   The  dis- 
cussions of  the  paper  demonstrated  the  debate  at  the  time 
about  the  causes  of  the  forces. 

The  study  of  interaction  forces  is  important  for  at 
least  four  reasons.   First,  they  contribute  to  collisions. 
The  understanding  of  these  forces  is  necessary  to  help  pre- 
vent future  occurrences.   Second  is  for  the  control  of 
ships  moving  in  proximity  to  each  other,  as  done  by  several 
Navys  to  accomplish  refueling  of  ships  at  sea.   The  third 
implication  is  in  the  design  of  future  ships  and  their 
operating  procedures.   The  adverse  effects  of  ship  inter- 
actions may  be  reduced  by  modifying  design  parameters  and 
effective  design  and  use  of  controls.   Finally,  for  those 
ships  in  operation,  the  present  rules  for  navigation  do 
not  recognize  the  interaction  effects.   Future  study  may 
show  the  need  to  modify  these  rules  so  as  to  provide  safer 
movement  of  the  merchant  fleet. 

With  this  perspective,  the  research  effort  here  at- 
tempts to  determine  the  extent  of  the  past  efforts  to  de- 
velop a  generalized  theory  for  ship  interactions  which 
would  be  practical  to  the  engineer.   The  view  is  toward  a 
general  model,  in  particular  one  which  could  be  used  as 
part  of  a  generalized  program  for  predicting  ship  trajecto- 
ries, 
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1. 2  Ship  Interaction 

1,2.1  General 

What  is  meant  by  interaction  forces?  Any  force  that  is 
only  observed  when  one  or  more  ships  (or  other  objects  such 
as  wall,  etc. )  are  present  can  be  classified  as  an  inter- 
action force.   The  motion  of  a  ship  through  the  water  mani- 
fests itself  by  waves  and  the  wake  formed  and  by  a  moving 
pressure  disturbance  in  the  water.   The  moving  pressure 
field  will  be  identified  as  the  potential  disturbance  where 
the  waves  will  be  referred  to  as  the  free  surface  disturb- 
ance. 

The  interactions  between  ships  operating  on  the  surface 
of  the  ocean  result  from  the  influence  of  both  the  free  sur- 
face and  potential  disturbances,  whereas  each  appears  to 
create  forces  by  different  mechanisms  and  differ  in  their 
spatial  distribution.   The  potential  field  generates  a 
distribution  of  pressures  in  the  fluid.   The  waves  also 
generate  a  pressure  field.   However,  this  rapidly  attenuates 
with  depth.   The  waves  also  change  the  water  surface  and 
this  introduces  forces.   If  the  frequency  is  low,  the  long 
duration  of  the  force  will  cause  significant  ship  motions. 
An  example  is  in  the  case  of  following  seas  where  direc- 
tional instability  may  occur  and  result  in  a  large  response 
of  the  ship. 

Very  close  to  the  ship,  the  potential  field  dominates 
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the  fluid  flow  pattern,  where  much  further  away  the  waves 
produce  the  flow  disturbances.   The  induced  flows  act  on 
other  bodies  (ships)  placed  in  the  perturbed  stream  to 
create  interaction  forces. 

1.2.1  Taylor  Results 

One  of  the  first  scientific  experiments  performed  and 
reported  on  the  interaction  between  ships  was  conducted  by 
David  v7.  Taylor  at  the  Model  Basin  in  V/ashington,  D.  C. 
The  results  are  reported  in  the  1909  Transaction  of 
Society  of  Naval  Architects  and  Ivlarine  Engineers,   A  good 
summary  of  the  work  can  be  found  in  Reference  2.   Here 
attention  is  focused  on  the  quantitative  measurements. 
The  results  are  presented  for  the  two  model  trials.   The 
first  is  with  relatively  thin  models  and  the  second  is  with 
fuller  models,   For  convenience,  the  model  characteristics 
are  also  presented.   In  each  case,  the  forces  are  reported 
for  only  one  of  the  models  as  indicated. 
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The  original  results  were  non-dimensional ized  by  the 
resistance  of  the  model.   Taylor  comments  that  the  results 
may  have  inaccuracies  due  to  problems  encountered  in  towing 
the  models  on  a  parallel  course,   He  also  stated  that  the 
results  for  the  thinner  models  were  less  consistant.   It  is 
difficult  to  say  if  the  data  was  erratic  or  if  the  results 
were  counter  intuitive,  thus  indicating  to  Taylor  that  they 
might  be  erratic.   Whatever  the  meaning,  both  results  are 
presented  and  evaluated  using  the  proposed  theoretical 
models.   Appendix  E  contains  calculations  of  the  possible 
errors  in  the  data  that  might  have  resulted  from  misalign- 
ment in  the  models. 

The  original  data  was  modified  to  develop  the  total 
forces  and  moments  on  the  models.   Also  a  change  in  the 
non-dimensionalization  was  made.   To  accomplish  these 
changes,  the  force  acting  on  the  model  was  determined  by 
multiplying  the  data  by  the  resistance  of  the  model.   Taylor 
states  that  the  models  were  towed  at  speed  between  tv/o  and 
three  knots.   A  value  of  2.5  knots  was  assumed.   The 
total  force  on  the  model  is  determined  by  adding  the  forces 
that  were  measured  at  each  end  of  the  model.   The  moments 
were  calculated  about  the  midship  position.   The  following 
curves  depict  Taylor's  data  after  the  indicated  changes. 

The  curves  are  drawn  for  the  ships  in  various  rela- 
tive positions.   The  abscissa  is  scaled  in  terms  of  the 
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FIGURE   2 

TAYLOR    MODEL  TESTS 


-3 


,3  xio 


a  X  lo 


1  X  lO 


-3 


-2  <  10 


—  :  -  +  xio 


i  X   CO 


-?.  xio" 


4  mo" 


-£,  ".  lo" 


1.0 


-.*  -.I- 

ASTERN 


,2  .'t-  .(, 

AHEAD 
SEPARATION      {^^O  15 


model  length.   The  attraction  and  repulsion  forces  are  self- 
explanatory  and  are  applicable  for  the  model  passing  on 
the  starboard  and  port  side.   The  moments  are  depicted 
as  being  either  bow  in  or  bow  out.   This  indicates  the  di- 
rection in  which  the  moment  is  acting  on  the  model.   The 
Thin  Ship  data  represents  Model  Number  83^  as  it  passes 
Model  Number  838.   The  Full  Ship  data  represents  Model 
Number  858  passing  Model  Number  866, 

1.2.3  Newton  Results 

Experiments  were  carried  out  at  the  Admiralty  Experi- 
mental Works  (AEV/),  Haslar,  England,  between  19^6  and  19^8 
to  test  out  the  feasibility  of  the  replenishment  of  war- 
ships with  fuel  tankers.  These  results  were  reported  at 
the  first  Symposium  on  Ship  Maneuverability  at  the  David 
Taylor  Model  Basin  in  Washington,  D.  C.  in  i960  by  R.  N. 
Newton,  The  tests  conducted  included  an  experiment  with 
constrained  models  of  the  Battleship  H.M.S.  KING  GEORGE 
and  the  R.F.A.  OLNA.  (Ship  A  and  B  Respectively) 

A.E.W.  used  freely  propelled  and  controlable  models  to 
study  the  behavior  during  maneuvering  alongside  and  breaking 
away.   The  work  also  contains  some  data  resulting  from 
tests  conducted  at  sea  with  ships.   The  results  are  presen- 
ted here  for  the  constrained  models.   The  force  moment  act- 
ing on  each  ship  has  been  converted  to  a  position  amidship. 
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FIGURE   3 


AE.W.     MODEL    TESTS 
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1 , 3  Previous  Analytic  Models 
1.3.1  Ellipsoid 

Sir  H.  T.  Havelock  proposed  a  model  of  ship  interac- 
tions based  upon  potential  flow.   As  an  example,  he  cited 
the  ellipsoid.   In  the  discussion  to  Reference  8,  the  solu- 
tions were  given  for  two  unequally  sized  ellipsoids  moving 
parallel  in  the  fluid.   The  forces  were  calculated  using 
Lagally's  theorem. 

The  following  singularity  distribution  describes  an 
ellipsoid  in  a  uniform  flow. 

SM(X)=2«C*Vhx 
WHERE 

1/C  =  ^  e/(l-e2)   -  2  log  [  ^1+  e)/  (i_e)  ] 

V  =  Speed  in  FT/SEC 

X  =  Distance  from  midbody  position 

e=  Ellipsoid  eccentricity 

Appendix  D  contains  sample  distributions  for  the 
ellipsoid  model.   Also,  the  results  for  the  sample  case 
cited  by  Havelock  are  recalculated  in  Appendix  F. 
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1.3.2  Rankine  Ovoid 

The  same  model  proposed  by  Havelock  for  use  with  the 
ellipsoid  could  be  applied  to  any  other  body  shape  for 
which  the  singularity  distributions  are  known.   The  Rankine 
Ovoid  is  one  of  the  simplest  bodies  to  generate  from  sin- 
gularities in  the  fluid.   It  is  formed  by  placing  two  equal 
but  opposite  singularities  in  a  uniform  flow  so  that  their 
axis  is  in  line  with  the  incoming  flow  and  the  source  is 
upstream  of  the  sink.   Given  the  singularities,  it  is  easy 
to  determine  the  size  of  the  body  that  is  represented  in 
the  fluid,   Silverstein  in  Reference  1^  developed  the  force 
on  a  Rankine  Ovoid.   Appendix  F  contains  the  recomputation 
of  a  case  demonstrated  in  his  paper.   Even  though  the  dis- 
tribution is  relatively  simple  and  the  body  shape  can  be 
determined,  it  takes  some  effort  to  solve  the  inverse  prob- 
lem.  The  desired  singularity  strengths   and  positions  for 
a  given  shape  requires  the  simultaneous  solution  of  two 
non-linear  equations. 

The  following  equations  describe  the  Rankine  Ovoid: 
^.(Strength)  (Dist)  (L/2) 


V  =  

L/2)2   /-^— ^2^2 

(2/2)2=  '' 


[(L/2)2  -  (Dist)2  ] 

k.    (Strength)  (Dist) 
V  J{B/2)^   +  (Dist) 2  ■ 

WHERE 

V  =  Speed  in  ft/sec. 

Dist  =  Distance  from  the  modbody  that  singulari- 
ties are  located. 

Strength  =  Singularity  Strength 
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The  solution  of  the  simultaneous  equations  was  accom- 
plished on  the  computer  .using  s  system  provided  program 
called  Zeroin,   A  copy  of  the  instruction  sheet  for  this 
program  is  included  in  Appendix  H.   A  tabulation  of  Ran- 
kine  Ovoid  distributions  is  provided  in  Appendix  D. 

The  experience  with  the  system  routine  that  solved  the 
simultaneous  equations  should  help  future  investigations. 
The  Zeroin  Program  is  one  of  two  that  were  available.   This 
particular  one  was  selected  mainly  because  it  offered  an 
automatic  feature  which  permitted  retrieval  of  intermediate 
steps  in  the  solution  so  that  the  program  results  could  be 
verified.   It  was  found  that  an  approximate  solution  to 
the  problem  using  the  semi-infinite  body  v/ould  provide  a 
good  starting  solution  for  the  program.   The  routine  would 
not  converge  to  the  problem  solution,  however,  until  the 
magnitude  of  the  initial  solution  vector  was  reduced  in 
size  by  two  orders  of  magnitude.   This  prevented  the  prog- 
ram from  jumping  past  the  solution  and  attempting  to  find 
another  solution  at  very  large  values  of  the  variables. 

1.3.3  Slender  Body  Approximation 

Newman  has  derived  the  forces  and  moments  acting  on  a 
slender  body  of  revolution  moving  parallel  to  a  wall.   In 
his  derivation,  it  was  assiimed  that  the  fluid  was  ideal  and 
of  infinite  extent,  thus  ignoring  any  free  surface  affects. 
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In  order  that  the  boundary  condition  at  the  wall  be  satis- 
file,  an  image  singularity  distribution  was  provided.   This 
is  analogous  to  two  bodies  of  revolution  moving  parallel 
and  abreast  through  an  infinite  ideal  fluid. 

The  approach  used  was  to  determine  the  singularity 
distribution  that  would  develop  fluid  boundaries  representing 
the  rigid  body  moving  through  the  mediiim.  The  singularities 
were  sized  using  a  two  dimensional  model  and  then  trans- 
formed to  three  dimensional  by  assuming  that  the  body  was 
slender  (rQ<<  L).   The  location  of  the  singularities  was 
determined  by  satisfying  the  linearized  boundary  condition 
on  the  body  surface.   It  was  pointed  out  tliat  for  the  body 
close  to   the  wall,  the  singularity  distribution  could  be 
approximated  by  a  general  bipolar  distribution.   This  re- 
sulted in  the  singularities  being  placed  along  the  length 
of  the  body  in  such  a  v/ay  as  to  be  displaced  from  the  center 
axis.   This  formulation  requires  the  wall  be  close 
(Z/L  «  1).  The  displacement  of  the  singularities  is  des- 
cribed by  the  following  formulas,  where  the  displacement 
from  the  centerline  axis  is  equal  to  (Zq  -  a). 

a^  =  Zo  ^-  ro^ 

Once  the  singularity  distributions  are  known,  the 
forces  acting  on  the  body  are  calculated  using  Lagally's 
theory.   It  is  easy  to  verify  that  the  singularity  distri- 
bution remains  within  the  body  and  thus  is  amenable  to  a 
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Lagally  force  calculation.   The  following  diagrams  show  the 
distributions  for  a  sample  ship  form. 

The  procedure  developed  the  singularity  distribution 
using  a  slender  body  assumption.   The  resulting  force  in- 
tegrals were  then  modified.   An  approximation  to  the  inte- 
gral is  based  upon  the  assimption  that  the  part  of  the  in- 
tegrand which  represents  the  singularity  as  a  function  of 
length  can  be  treated  as  a  constant  value. 

Nevmian's  v/ork  produced  the  following  results  for  the 
problem. 
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WHEPffi   Ze  =  Separation  of  the  body  axis  from  the  wall 
To  =  Radial  dimension  at  position  X  along  the 
body 

Newman  also  presented  a  sample  case  to  demonstrate  the 
results.   The  body  was  assumed  to  be  an  ellipsoid  des- 
cribed in  the  following  form: 

-2 


\mERE 


r,  (x)  =  r^  *(1-^  5^7  2) 


r^  =  Maximum  radius  of  the  body 
X  =  Longitudinal  position  measured  from  the 
midship  station. 
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Substituting  into  the  Force  Equation, 

WHERE 

Q^  is  the  Legender  Function  of  the  Second  Kind 

AS   ^  -  I 

'  -2  MAX  ' 

The  force  integrals  were  evaluated  for  this  case  and 
the  results  presented  as  follows. 


AND    ^  =  ^*4^ 
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In  the  concluding  comments  of  the  paper,  the  author 
states,  " , . . . it  seems  likely  that  for  bodies  with  full  or 
blunt  ends  this  [moment]  integral  will  be  dominated  by  the 
large  values  of  S'(x)  at  the  two  ends,  and  thus  the  pitch 
moment  will  be  positive  (bow  down)  if  the  bow  is  blunt 
compared  to  the  stern,  and  vice  versa".   In  this  comment, 
the  pitch  is  analogous  to  the  yaw  in  the  previous  defined 
orientations.   He  is  predicting  a  'bov;  out'  condition  if 
the  sterns  are  blunt  compared  to  the  bows.   The  author 
also  notes  that  the  force  is  always  attractive  and  ".,,  for 
geometrically  affine  bodies,  having  the  same  cross  sec- 
tional distributions  but  different  lengths,  the  force  will 
be  inversely  proportional  to  the  length  whereas  the  moment 
will  be  independent  of  the  length",   From  this  result, 
the  nondimensional  value  of  maximum  force  is  inversely 
proportional  to  the  cube  of  the  L/B  ratio. 
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CHAPTER  II. 
FORCE  MO  DEI. 

2.1  Model  Descriiotion 

The  proposed  model  for  predicting  ship  interaction 
forces  uses  the  hypothesis  proposed  by  Havelock  that  the 
potential  field  could  explain  the  ship  interaction  pheno- 
mena.  This  approach  seems  reasonable  for  several  reasons. 
First,  the  ship  moving  at  the  water's  surface  produces  a 
flow  field  similar  to  that  produced  in  an  infinite  fluid. 
Second,  the  potential  field  dominates  the  flow  near  the 
body.   Pressures  in  this  region  fluctuate  spatially  and  the 
changes  find  relief  on  the  surface.   The  movement  of  this 
pressure  field  causes  the  wave  system.   Third,  for  most 
cases  of  interest,  the  ships  are  operating  below  their  de- 
sign speed.   This  is  particularly  true  for  naval  ships  re- 
plenishing at  sea.   The  fine  entry  of  these  ships  causes 
minimal  wave  disturbances  when  at  replenishment  speeds. 

The  potential  flows  without  a  free  surface  can  be 
modeled  using  bodies  moving  in  an  infinite  fluid  with  image 
systems  to  maintain  the  flow  conditions  on  the  waterline 
plane.   This  condition  requires  that  there  be  no  normal 
component  of  flow  at  the  interface  representing  the  undis- 
turbed water  level.   A  body  of  revolution  moving  horizon- 
tally in  the  fluid  has  a  symmetry  of  flow  and  is  used  to 

model  the  ships.   The  maximum  diameter  of  this  body  is  taken 
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as  the  ship's  beam  because  of  the  major  influence  this  di- 
mension has. 

The  force  acting  on  the  total  body  of  revolution  would 
represent  twice  the  force  experienced  by  the  ship.   The 
force  results  from  the  integration  of  pressures  on  (the 
lower)  half  of  the  body  of  revolution.   The  symmetry  of 
flow  will  insure  the  same  pressures  and  thus  forces  on  the 
other  half  body. 

The  major  reason  for  selecting  a  body  of  revolution  was 
the  ability  of  using  an  axial  distribution  of  singularities 
to  describe  the  flov;.   Hess  et,  al,  have  shown  the  compu- 
tations required  for  representing  an  arbitrary  body  in  an 
ideal  fluid  by  a  surface  distribution  of  singularities. 
The  computational  effort  for  each  ship  would  be  much  greater 
for  surface  distribution  and  as  such  is  not  desirable  for 
use  in  a  ship  trajectory  model. 

In  total,  four  models  were  developed  and  programmed. 
The  first  two  incorporated  the  ellipsoid  and  Rankine  Ovoid 
as  discussed  earlier.   The  third  and  final  one  involved  a 
modified  ellipsoid  and  slender  body  respectively,  each  with 
cross  flow  corrections.   The  first  two  used  the  derivations 
of  Havelock  and  Silverstien.   The  third  is  only  a  special 
case  of  the  fourth  model.   The  slender  body  model  will  be 
discussed  in  more  detail.   Appendix  I  contains  a  logic  flow 
chart  of  the  computer  program  for  this  model, 
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The  slender  body  model  uses  bodies  of  revolution  to 
model  the  ships  as  described  previously,  the  axial  distri- 
bution of  sources  and  sinks  is  obtained  by  assuming  the 
body  of  revolution  is  slender.   In  two  dimensions,  the  sin- 
gularities are  sized  using  the  follov/ing  equations.   The 
derivation  for  three  dimensions  are  presented  in  Appendix  C, 
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This  axial  distribution  of  singularities  determines  the 
body  in  a  uniform  flow.   When  another  body  comes  near,  the 
flow  field  is  changed,  so  that  it  is  no  longer  uniform. 
As  the  bodies  approach  each  other,  their  affect  in  altering 
the  flow  field  causes  a  deformation  in  the  fluid  boundary 
which  describes  the  rigid  bodies.   To  handle  this  problem 
an  additional  distribution  of  dipole  singularities  are 
placed  along  the  axis.   Thes   ire  oriented  in  the  plus  or 
minus  direction  on  the  body  axis  and  are  sized  to  counter 
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the  effects  of  local  cross  flow  on  the  bodies.   The  change 
in  the  flow  in  the  axial  direction  is  handled  by  perturbing 
the  original  singularity  distribution. 

The  boundary  condition  on  each  body  gives  no  flov/ 
across  the  boundary  which  describes  the  rigid  body.   In  our 
case,  this  is  a  body  of  revolution.   At  each  location  along 
the  body  there  is  a  radial  distance  at  which  the  velocity 
(radial)  must  be  zero.   If  the  cross  flows  were  uniform, 
the  motion  of  an  infinite  cylinder  could  be  represented 
exactly  in  the  fluid  by  placing  the  appropriate  doublet 
distribution  on  the  axis.   The  flow  pattern  around  a  slender 
body  of  revolution  can  be  approximated  locally  by  an  infi- 
nite cylinder  because  the  change  in  diameter  would  be  small 
for  a  change  in  the  axial  direction.   For  the  distances 
between  bodies  of  interest,  it  was  estimated  that  the  uni- 
form flow  assumption  for  any  one  segment  of  the  body  would 
be  nearly  correct.   The  dipoles  were  sized  with  flov/  of  the 
value  experienced  on  the  body  axis  as  a  measure  of  the 
uniform  flow  strength.   For  the  actual  flow  conditions,  the 
boundary  location  is  not  calculated.   It  is  assumed  that 
the  deviation  from  the  exact  position  would  be  small.   The 
dipole  distribution  is  sized  by  the  following  rule.    See 
Appendix  C  for  derivation. 


50   iy)     =     -^  '  V(')  '  ^^^^^^  / ^^ 
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If  the  perturbed  flov/  is  known,  the  changes  in  the 
singularity  distribution  can  be  readily  determined.   For 
the  case  of  two  bodies  near  each  other,  the  change  in  the 
singularity  distribution  in  one  of  the  bodies  is  reflected 
in  a  slight  change  in  the  perturbed  flow  field  experienced 
by  the  others.   V/hat  is  desired  is  a  simultaneous  solution 
to  the  problem.   This  is  achieved  by  performing  several 
iterations  of  the  flow.   Taking  first  one,  than  the  other 
body,  change  its  singularity  distribution  to  counteract 
the  perturbed  flow,  and  then  determine  the  resulting 
affect  on  the  flow  field  at  the  other  body.   Changes  in 
the  perturbed  flow  were  insignificant  after  about  three  to 
five  iterations  except  when  the  bodies  were  etremely  close 
and  then  it  took  only  about  ten  iterations  to  achieve  con- 
vergence to  the  solution. 

2.2  Interaction  Force  Calculations 

The  interaction  forces  acting  on  a  rigid  body  moving 
in  an  ideal  fluid  v/ithout  viscosity  can  be  determined  if 
the  fluid  velocity  field  can  be  represented  by  a  potential 
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function.   In  principle,  the  method  for  determining  a  force 
on  any  body  would  involve  the  calculation  of  the  pressure 
on  the  body  surface  followed  by  the  integration  to  obtain 
the  total  force.   This  was  first  accomplished  by  either 
Munk  or  Lagally  for  bodies  generated  with  singularities. 
Lagally  noted  that  for  a  body  placed  in  a  steady  stream 
this  integration  took  a  form  that  would  be  analogous  to  a 
derivation  which  treats  a  force  as  acting  between  the  sin- 
gularities.  The  apparent  force  between  a  pair  of  singulari- 
ties is  determined  in  three  dimensions  by  the  following 
formula. 


f:^  =  ^TT^m.^yii  /  2 


Where  a  positive  force  indicates 
an  attraction. 

The  total  force  would  be  determined  by  adding  all 
pair  contributions.   The  force  acting  on  a  body  which  is 
represented  by  a  continuous  source  distribution  along  its 
axis  has  the  follov/ing  integral  form. 

Where  [a.b]    and  [a',b']   indicate 
the  line  segments  with  singularities. 

The  results  of  Lagally  were  expanded  by  Ciimmins  to 
include  the  forces  between  doublet  singularities  for  the 
general  time  dependent  case.   His  derivation  involved  the 
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following  assumptions. 

1,  The  velocity  field  is  irrotational 
and  has  a  velocity  potential  ?  (x,y,z,t) 

2,  If  the  body  were  not  present,  the 
stream  would  have  a  velocity  potential  j^, 
which  we  call  the  potential  of  the  "undis- 
turbed stream". 

3,  There  are  no  singularities  of  the 
undisturbed  stream  in  the  region  occupied 
by  the  body. 

^■.      The  boundary  condition  at  the  sur- 
face of  the  body  is  satisfied  by  super- 
imposing a  system  of  singularities  upon 
the  undisturbed  stream,  such  singularities 
falling  v;ithin  the  region  which  the  body 
would  occupy.   The  potential  of  the  system 
of  singularities  is  designated  by  ^•^.      Then 

J  =  ^  +  ^b 

Cummins  found  that  the  forces  and  moments  calculated 
were  of  the  same  form  as  expressed  in  the  Lagally  deriva- 
tion, but  with  additional  terms.   He  expressed  the  Lagally 
force  term  as  a  product  between  local  singularities  and 
the  total  flow  perturbation  rather  than  as  the  sum  of  in- 
dividual contributions  from  each  singularity.   The  two 
representations  for  the  force  are  identical.   Zucker  gives 
a  svimmary  of  the  results  and  indicates  the  contribution 
by  Landweber  and  Yih  in  substantiating  Cummins'  results 
and  evaluating  the  moment  term  resulting  from  the  change 
in  the  fluid  flow  with  time. 

The  force  on  the  body  will  be  as  follows  for  Cummins' 

assimptions. 
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^  St 

WHERE 

¥"  ^^  Volume  of  the  body 
Vo  -  Is  the  translational  velocity  of  the  body 

v;  -    Is  the  location  of  the  centroid  of  the 
''    body  relative  to  the  origin  of  the  moving 
system. 

s '  The  force  if  the  body  were  not  rotating 
and  the  undisturbed  stream  were  steady 
(lagally  ?orce) 

F. -Force  due  to  change  of  the  flow  with  time 


Fj-Force  arising  from  body  rotations 


M,  -Moment  if  the  dlow  were  steady  (Lagally) 


1-^  -  Moment  due  to  the  flow  changing  with 
time 

h^  '  Moment  effect  due  to  rotation  of  the 
body 

FOR  SITUATIONS  DESCRIBED  BY  ONLY  SOURCES  AND  DOUBLETS, 
F^  AND  M^   ARE  ZERO. 
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FOR  A  SOURCE  OF  STRENGTH 


F,-  -^fTf  -^^ 


WHERE    ^'   defines  the  fluid  velocity  in  the 
undisturbed  stream 


J-      V 


FOR  A  DOUBLET  0^  MOMENT 


NOT  EVALUATED  BY  CUMMINS  ARE  THE  MOMENT  DUE  TO  THE 


FLOW  CHANGING  WITH  TIME. 
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^'  'H  f(.'f)  iCr^^)  J(r 


This  term  was  evaluated  by  Landwebber  and  Yih, 
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The  force  calculations  programraed  on  the  computer 
used  at  least  one  of  the  forms  discussed  previously.   As 
more  complicated  and  detailed  models  were  developed  the 
form  of  the  force  calculation  was  selected  to  improve  com- 
putational efficiency.   For  the  Rankine  Ovoid  and  ellipsoid 
which  do  not  contain  cross  flow  corrections,  the  simplest 
statement  of  Lagally's  force  was  used.   When  it  became 
necessary  to  introduce  dipoles  to  correct  for  the  cross 
flow,  the  final  form  of  Cummins  was  used.   The  final 
slender  body  model  takes  into  account  the  terms  that  are 
represented  by  Fj^  and  Mj_ .    This  represents  the  non-zero 
terms  of  the  interaction  force  for  steady  parallel  motion 
of  the  bodies. 

The  above  equations  relate  to  a  three  dimensional 
flow  problem.   The  two  dimensional  case  is  similar  with 
the  Lagally  force  being  of  the  following  form: 

Havelock  developed  in  Reference  7  the  force  and  moments 
acting  on  a  two  dimensional  body  by  using  contour  integra- 
tion. 

FOR  SOURCE  DISTRIBUTION 
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WHERE 

1.  Suffix  "S"  refers  to  the  given  distribution 

in  the  liquid,  and  "r"  to  the  image  system 
within  the  surface  of  the  body. 

2.  Summation  extends  over  the  external  and  inter- 

nal sources  taken  in  pairs. 


FOR  DIPOLE  DISTRIBUTION 


Where:     —  the  doublet  is  of  moment  Nl  and  makes 
an  angle  c{       with  the   x   axis  (ship). 

And;   ^  being  the  angle  between  K°    and  the  vector 
drawn  from  the  doublet  (r)  to  the  doublet  (s). 

The  total  moments  results  from  moments  of  forces  already 
determined  plus  an  additional  term. 
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2. 3  Computation 

The  earlier  discussion  of  the  slender  body  model  gave 
the  general  steps  involved  for  determining  the  interaction 
forces.   In  this  procedure  there  were  two  major  areas 
requiring  considerable  effort.   First,  it  was  necessary  that 
derivatives  of  the  flow  potentials  be  derived  for  the  force 
calculations.   These  potentials  and  their  derivatives  were 
defined  in  terms  of  a  moving  axis  system  of  the  ship  in 
which  the  forces  were  acting.   The  necessary  transformations 
and  derivations  are  developed  in  Appendix  B, 

The  second  area  involved  the  determination  of  the  total 
force  acting  on  the  ship.   The  slender  body  was  defined  by 
a  continuous  distribution  of  singularities.   The  resulting 
flows  in  the  fluid  v/ere  determined  by  integrating  the  con- 
tributions from  the  singularity  distribution.   Also  the 
forces  were  determined  by  integrating  the  individual  force 
contributions.   In  both  cases,  the  numerical  integration 
used  values  of  the  integrand  at  the  same  location  on  the 
body  axis.   These  coincided  with  the  stations  at  which 
the  sectional  area  curve  had  been  defined.   The  integration 
was  performed  in  each  case  using  Simpson's  Rule  with  ai) 
odd  number  of  stations  (21),  evenly  spaced  along  the  body 
axis. 
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CHAPTER   III. 

COr/JFARISON   OF    RESULTS 

3.1   Taylor  Model  Tests 

The  Rankine  Ovoid  and  ellipsoid  models  were  developed 
and  checked  with  the  available  results  then  used  to  model 
the  forces  and  moments  reported  by  Taylor.   Figure  1    de- 
picts the  results.   These  predictions  differ  in  two  ways 
with  Taylor.   First,  they  do  not  predict  the  secondary  re- 
versal in  the  moment  that  is  evident  at  longitudinal  se- 
parations of  .6*1.   Secondly,  the  largest  moment  is  pre- 
dicted for  longitudinal  separations  above  1.0  xL.  This 
variation  can  be  attributed  to  the  deformation  of  the  boun- 
daries of  the  bodies  as  they  approach  each  otlna?.   If  the 
boundary  was  maintained  in  the  original  position,  the  flow 
would  have  been  restricted,  thus  increasing  the  fluid 
velocity  and  contributing  to  a  lower  pressure  on  the  bodies. 
This  reduction  in  pressure  would  contribute  a  moment  that 
would  cause  a  reversal  in  sign.   The  inability  of  the  pre- 
vious theoretical  models  to  explain  the  reversal  in  sign 
of  the  moment  curve,  led  to  the  development  of  a  model  which 
incorporated  cross  flow  corrections.   The  correction  for 
cross  flow  indicates  the  reversal  in  the  moment  can  be 
modeled.   However,  the  magnitudes  were  not  in  full  agree- 
ment with  those  obtained  by  Taylor.   Computations  were  made 
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at  higher  speeds  with  no  variation  in  the  results.   This 
was  expected,  th^   plot  involved  the  non-dimensional  moment 
coefficient.   For  closer  distances,  the  predicted  ampli- 
tude of  the  moments  increased  only  slightly.  (The  model 
tests  were  conducted  at  a  very  small  separation. ) 

The  dimensions  of  the  models  used  in  Taylor's  experi- 
ments were  inputs  to  the  slender  body  model.    Calculations 
were  made  for  the  same  cases  as  described  by  Taylor.   The 
following  curves  show  the  results  plotted  for  comparison 
purposes. 

Additional  curves  are  drawn  which  show  the  forces  pre- 
dicted for  the  Taylor  full  bodied  models,  as  the  lateral 
separation  is  increased.   The  reversal  in  moment  diminishes 
as  the  ships  separate  because  there  is  no  longer  a  restric- 
tion in  the  flov/  to  cause  a  pressure  drop  between  the  ships. 
Beyond  a  distance  of  2  beams, the  results  agreed  with  those 
obtained  v/ith  a  model  which  did  not  incorporate  cross  flow 
corrections.   Only  beyond  a  separation  of   4  beam  widths 
interaction  effects  were  negligable. 

The  model  tests  conducted  by  Taylor,  and  the  slender 
body  predictions  differ  in  several  ways.   Some  of  the  more 
significant  assumptions  used  for  the  slender  body  were: 

1,  The  potential  field  was  assumed  most 
important  (No  free  surface) 

2.  The  ships  were  represented  by  a  body 
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FIGURE   9 
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of  revolution,   (Lessen  effective  beam) 

3.   Squat  and  sinkage  were  not  incorpo- 
rated. 
Other  differences  in  the  results  may  be  attributed  to  either 
asemmetry  in  the  models  or  inaccuracies  in  the  measurement 
techniques.   The  apparatus  used  by  Taylor  was  not  indicated, 
V/hen  the  tests  were  conducted,  the  resistance  measurements 
used  a  system  of  weights,  pulleys  and  levers.   Taylor  even 
commented  on  the  problem  of  conducting  the  model  tests. 
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3. 2  Comparison  of  Results  with  Newman 

In  order  that  a  comparison  be  made  between  the  model 
developed  and  the  results  of  Newman.   The  sample  problem 
in  Newman's  paper  was  used.   The  sectional  area  curve  for 
the  ships  moving  in  parallel  were  determined  using  the 
radius  formula  presented  earlier.   Two  ships  moving  in 
parallel  were  defined  with  the  same  geometry.   The  initial 
calculation  assumed  a  length  of  600  feet,  a  beam  of  100 
feet  and  a  speed  of  10  knots.   The  lateral  displacement 
was  determined  and  plotted  according  to  the  non-dimensiona- 
lization  used  by  Newman. 

Note  that  Zq/v^  is  equivalent  to  the  separation  dis- 
tance between  the  axis  of  the  two  bodies  divided  by  their 
beam  dimension.  Other  beam  to  to  length  ratios  were  cal- 
culated by  varying  the  length  of  the  bodies  while  main- 
taining the  beams  constant.  These  included  the  L/B  of  6, 
10,  20,  6o,  The  total  force  on  the  body  of  revolution 
determined  using  the  model  was  non-dimensionalized  by  the 
maximum  value  determined  by  Nev/man.  The  following  graphs 
depict  the  results  developed  in  this  thesis. 

The  first  graph  shows  the  results  for  a  L/B  of  6  for 
the  various  theoretical  models.   Calculations  were  then 
made  for  the  resulting  forces  with  and  without  the  cross- 
flow  correction  which  used  the  dipole  distribution.   It 
was  expected  that  the  results  of  Newman  and  the  proposed 
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model  would  approach  each  other  in  the  limit  of  high  L/B 
ratios  and  for  very  slender  bodies  because  the  singularity- 
distributions  would  be  nearly  identical  and  lie  very  close 
to  that  used  by  Newman. 

The  next  graph  shows  the  change  in  the  non-dimensiona- 
lized  maximum  force  as  a  function  of  the  length  to  beam 
ratio.   Note  that  the  singularity  model  approaches  the  New- 
man result  in  the  limit.   On  the  graph,  the  error  of  the 
computation  is  indicated.   The  solution  with  the  dipole 
correction  seems  to  converge  initially  and  then  for  very 
high  ratios  diverge.   This  tends  to  indicate  that  the  di- 
pole correction  may  be  useful  for  all  ratios  of  interest. 

The  error  introduced  during  integration  was  easy  to 
approximate  for  the  axial  distribution  of  singularities 
because  the  fourth  derivative  of  the  integrand  should  be 
exactly  zero.   However,  for  the  case  with  the  dipole  cross 
flov/  correction  this  may  not  be  the  case  and  it  is  expected 
that  the  introduced  error  would  be  larger.   This  may  ex- 
plain the  divergence  of  the  solution  for  high  L/B  (large 
lengths).   See  Appendix  E.  for  derivation  of  error  estimate. 

Another  graph  shows  more  directly  the  variation  in 
the  results  for  the  various  length  to  beam  ratios.   In  this 
graph  the  shaded  area  indicates  the  effect  of  the  dipole 
correction. 

The  last  graph  is  a  compilation  from  various  sources 

of  interaction  forces.   They  are  plotted  in  the  same 
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non-dimensional  terms.   Tor  the  ships,  the  forces  were 
doubled  to  approximate  the  forces  that  would  be  observed 
between  bodies  of  revolution.   The  non-dimensionalization 
with  respect  to  the  maximum  force  predicted  by  Newman  is 
consistant  with  the  other  graphs.   Plotted  also  are  the 
solutions  predicted  for  a  body  of  revolution  with  L/B 
ratios  of  6  and  10  with  the  geometry  used  by  Newman.   The 
plot  indicates  that  the  model  tests  agree  more  closely  to 
the  results  predicted  by  the  proposed  model  with  the  results 
of  the  tests  by  Newton  showing  the  greatest  variation.   If 
the  results  had  been  mislabeled  the  results  would  have  been 
more  consistant  v/ith  the  proposed  model  predictions. 

The  basic  differences  in  the  proposed  model  and  that 
of  Newman  are,  first,  the  surface  condition  on  the  body  is 
accounted  for  by  displacing  the  singularities  from  the 
centerline  in  the  case  of  Newman,  v/here  the  proposed  model 
uses  a  dipole  distribution.   It  is  obvious  that  the  pro- 
posed model  will  satisfy  the  surface  condition  better  when 
the  induced  velocity  is  fairly  uniform  over  the  body  pro- 
file.  From  the  data  presented,  the  relatively  small  L/B 
and  the  separation  distances  of  interest  justify  the  use 
of  dipole  corrections.   A  second  difference  in  the  models 
is  in  the  evaluation  of  the  integral  solutions.   Newman 
uses  an  approximation  to  the  solution  which   emphasizes 
the  need  to  restrict  the  solution  to  the  slender  bodies, 

whereas   this    model  uses  niAmerical  integration. 
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CHAPTER  IV 
CONCLUSIONS 

The  major  effort  of  the  thesis  has  been  to  provide  a 
theoretical  model  for  the  prediction  of  forces  "between 
ships.   A  computer  program  has  been  developed  which  carries 
out  the  computations  for  ships  in  rectilinear  motion  using 
the  potential  flow  model  described.   The  ellipsoid  model 
proposed  by  Havelock  was  unable  to  predict  the  form  of  the 
results  as  indicated  by  empirical  model  tests.   These 
tests  demonstrated  a  reversal  in  the  moment  force  as  the 
ship  separation  decreased.   This  led  to  the  development  of 
the  slender  body  model  with  cross  flow  corrections.   The 
programs  which  were  developed  for  ships  in  rectilinear  mo- 
tion ignore  time  dependent  forces.   The  data  presented 
here  represents  the  interactions  of  ships  in  deep  water, 
on  parallel  courses  and  with  each  ship  at  the  same  speed. 
Orientation,  other  than  parallel,  would  have  resulted  in 
a  change  in  the  relative  positions  of  the  ships  with  time 
and  would  involve  other  force  terms  as  derived  by  Ciimmins, 
but  were  not  programmed  here. 

The  investigation  of  past  efforts  indicate  a  lack  of 
good  model  or  full  scale  data  to  test  a  theoretical  result. 
The  differences  experienced  in  the  analysis  of  Taylor's 
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results  could  not  be  fully  explained  without  more  empirical 
data.   This  can  only  be  obtained  with  ship  models  because 
of  the  parallel  course  requirements. 

In  meeting  the  objectives  of  the  thesis,  the  axial 
distribution  of  singularities  and  the  force  calculation 
techniques  provide  additional  inputs  to  a  general  trajec- 
tory model.   When  the  ship  approaches  other  bodies  or 
boundaries  in  the  fluid,  this  work  will  help  to  predict 
the  interactive  forces  so  that  more  inclusive  trajectory 
calculations  can  be  made.   This  is  applicable  to  the  deter- 
mination of  the  factors  involved  in  ship  collisions,  the 
prediction  o^  the  control  parameters  for  at  sea  replenish- 
ments, or  for  future  ship  design  characteristics. 

The  performance  of  model  tests  may  indicate  results 
that  differ  from.  Taylor's.   The  design  of  model  experi- 
ments must  take  into  account  the  large  forces  introduced 
by  any  misalignment  of  the  models.   Also,  the  wall  effects 
may  introduce  erroneous  readings.   The  computer  program 
has  been  written  to  include  more  than  two  ships  in  proxi- 
mity.  It  could  be  used  to  estimate  the  wall  effects  in 
the  tow  tank  experiments  by  introducing  the  first  few  image 
distributions  to  approximate  the  wall  boundary  conditions. 
The  following  figure  shows  the  orientation.   It  will  be 
important  that  some  tests  be  repeated  with  the  models  in 
a  reverse  geometry  to  insure  accurate  readings.   This  would 
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FIGURE   1^ 

THEORETICAL       CALCULATION      OF     WALL      EFFECTS 
ON        INTERACTION      FORCE    MODEL      EXPERIMENTS 
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give  a  measure  of  the  wall  effects  and  asymmetries  in  the 
models. 

The  program  can  be  extended  to  cover  shallow  water 
effects.   For  very  shallow  water,  a  tv/o  dimensional  model 
could  "be  developed  by  only  a  few  changes  in  the  present 
program.   These  would  involve  the  change  in  sizing  the 
singularity  distributions  and  in  the  Lagally  force  calcu- 
lation,  ^or  other  shallow  water  conditions,  the  boundary 
conditions  could  be  satisfied  by  a  series  of  image  distri- 
butions as  shown  in  Figure  15. 
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FIGURE  15 

SHALLOW     WATER      APPROXIMATIOM 
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CHAPTER  V. 

RECOMvTENDATIONS 

The  work  accomplished  in  the  thesis  points  out  the 
need  for  continued  study  of  the  interaction  betv/een  ships. 
The  following  list  indicates  possible  areas  for  future  in- 
vestigation. 

1.  Perform  model  experiments  to  obtain  data  to 
verify  Taylor's  results,  including  shallow  water  effects 
as  described. 

2.  Include  remaining  force  terms  in  order  to  calculate 
for  unsteady  ship  motions. 

3.  Incorporate  program  in  trajectory  model. 

k.      Determine  stable  positions  alongside,  and  test 
with  operating  experience. 

5.  Simulate  refueling  maneuvers  and  compare  with 
operator  experience. 

6.  As  a  final  test,  attempt  to  reconstruct  published 
collisions  involving  ship  interactions. 

7.  Investigate  feasibility  of  extending  the  theoreti- 
cal model  to  include  free  surface  condition  by  distributing 
singularities  as  used  by  Havelock  and  developing  the  re- 
quired theorems  to  evaluate  the  interaction  forces  as  accom- 
plished by  R.  D,  Zucker  for  vortex  sheets. 
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APPENDIX  A 


SSAC  (I,  J)  = 


Nondimensional  sectional  area  at  Station 
J  of  Ship  I 


SAC   (I, J) 
DSAC  (I, J) 

DISPL(I) 
LENGT(I) 
BEMl  (I) 
VELOC  (I) 
CP  (I) 
CM  (I) 
DIST  (I) 

S  (I, J) 

SM  (I, J) 

SD  (I, J) 

NN  (I) 

KK  (I) 

X  (I) 
Y  (I) 
CI  (I) 


=  Sectional  area  at  Station  J  of  Ship  I 

=  Derivative  of  sectional  area  at  Station  J 
of  Ship  I 

=  Displacement  of  Ship  I   (Not  used) 

=  Length  of  Ship  I 

=  Beam  of  Ship  I 

=  Velocity  (retilinear)  of  Ship  I 

=  Prismatic  coefficient  of  Ship  I 

=  Midship  coefficent  of  Ship  I 

=  Distance  from  midship  where  R,  0.  singu- 
larities are  placed  to  generate  Ship  I 

=  Initial  (no  interaction)  source  strength 
at  position  J  on  the  axis  of  Ship  I 

=  Source  strength  at  position  J  on  the  axis 
of  Ship  I 

=  Dipole  strength  at  position  J  on  the  axis 
of  Ship  I 

=  Number  of  stations  used  to  input  non- 
dimensional  sectional  area  curve  and  to 
perform  calculations 

=  Index  for  Ship  I  to  indicate  model  used 
to  calculate  forces  (not  used) 

=  X  position  of  Ship  I  in  the  fluid 

=  Y  position  of  Ship  I  in  the  fluid 

=  Heading  angle  of  Ship  I  in  the  fluid 

61 


PI  =  The  constant 

RON  =  Fluid  density 

M  =  Index   (not  used) 

Wl  =  Index   (not  used) 

N      =  Number  of  ships  interacting 

(Program  can  handle  up  to  5  with  space  allo- 
cation) 

XFl(J)  =  Component  of  force  in  the  X  direction  at 
Station  J 

YFl(J)  =  Component  of  force  in  the  Y  direction  at 
Station  J 

R^l(J)  =  Moment  at  Station  J  on  Ship  I 

XF2  =  Component  of  force  in  the  X  direction  of  Ship  2 

J  =  Index  for  indicating  position  along  Ship  I 

L  =  Index  for  indicating  position  along  Ship  K 

XI (J)  =  Longitudinal  distance  to  Station  J  on  Ship  I 

Y1(L)  =  Longitudinal  distance  to  Station  L  on  Ship  K 

I      =  Index  to  indicate  first  ship  of  the  pair  for 
force  calculation 

K      =  Index  to  indicate  second  ship  of  the  pair  for 
force  calculation 

AFX(I)  =  Total  force  on  Ship  (I)  in  X  direction  (Out- 
put on  non-dimensionalized  ED-form) 

AFY(I)  =  Total  force  on  Ship  (I)  in  Y  direction 

AM(I)   =  Total  moment  on  Ship  I 

ACM(I)  =  Total  non-dimensionalized  moment  for  Ship  I 
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QR   (I, J)   =  Total  flow  at  Position  J  of  Ship  I  (not 
used) 

QA   (I, J)   =  Total  axial  flow  at  Position  J  of  Ship  I 

QN   (I,  J)   =  Total  transverse  flov/  at  Position  J  of 
Ship  I 

DQA  (I, J)   =  Total  axial  contribution  of  flow  gradient 

at  location  J  of  Ship  I 
DQN  (I, J)   =  Total  transverse  contribution  of  flow 

gradient  at  Location  J  of  Ship  I 

QAZ  (L)    =  Axial  flov/  contribution  at  J  on  Ship  I 
due  to  singularity  at  L  on  Ship  K 

QNZ  (L)  =  Transverse  flow  contribution 

DGNZ  (J)  =  Transverse  gradient  contribution 

QRZ  (J)  =  Flow  contribution 

DQAZ  (J)  =  Axial  gradient  contribution 

DXI   (J)   =  Station  spacing  on  Ship  I  between  Station 

J  and  J+1 

DXK   (L)   =  Station  spacing  on  Ship  K  between  I  and 

L+1 

AEX  (I)   =  Focal  distance  of  ellipse  measured  from 

midship  position 
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APPENDIX  B 

Derivation  of  Axis  System  and 
Transformations  and  Potential  derivatives 
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DEFINITION    OF    AXIS    SYSTEMS 
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CONVENTION     FOR     DETERMING    -RELATIVE 

DISTANCES        BETWEEN         SINGULARITIES       IN 
THE      FORCE         CALCULATIONS     OF    SHIP     1 


t 

+Y 


SHI?    2 
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FORCE    ON     SHIP    i     DUE 
TO    SINGULARITIES       M,  AND  M, 


■*-  +x 


SHIP  2 


Yi*swlcia')] 


SHIP  1 


-Xi'  SwLcx  C-i>i 


ACi)  =  KX  -  XI'  cos  Cci  C-iTl  +  Yi  *  Cos  Cci  C2  )1 


F  =  (4Tr,f^m,hn^/p2)f 
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RESOLUTION      OF     FORCE     ON      SHIP  i 
INTO      AXIAL     AND      NORMAL       COMPONENTS 
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APPENDIX  D 


SINGUT.ARITY  DISTRIBUTIONS 


The  following  graphs  depict  some  of  the  representative 
singularity  distributions. 

Model  of  an  ellipse  with  I/B=6,  L=  600  ft.,  Y=  10  kts. 
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The  modi-Tied  ellipsoid  has  three  singularity  distri- 
when  placed  moving  parallel  to  a  wall  with 
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The  distributions  used  to  model  the  Taylor  full  ship 
results  are  as  follows.   The  situation  when  the  models 
are  abeam  with  a  3.9  ft.  separation  is  represented  (end 
corrections  included). 
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APPENDIX   F 

Presentation  of  Ellipsoid  and 
Rankine  Ovoid  Results  (Compare  with  original  work) 

The  results  are  the  same  except  for 
a  scale  change.   Havelock's  scale  should  be  divi- 
ded by  10,  where  Silverstien's  should  be  divided 
by  (^4Tr)\   or  158. 


88 


APPENDIX  G 
COMPUTATION  E^T^ICIE^CY 

The  programs  written  without  cross  flow  corrections 
take  only  a  fraction  of  a  second  to  evaluate  interaction 
forces  for  each  individual  ship  position.   The  slender 
body  model  with  cross  flow  corrections  takes  various  times 
for  execution  depending  on  the  number  of  iterations  of 
the  flows  required.   For  the  samples  rvin,  two  ships  with 
twenty  one  stations  each,  took  three  seconds  for  five 
iterations  increased  to  ten,  the  time  doubled  to  six  sec- 
onds for  each  ship  position  calculated. 
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Applications  Program  Series 
AP-21  revision  2 


APPENDIX  H 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
INFORT'IATION  PROCESSING  CENTER 


NAME: 


DESCRIPTION ; 


September  28,  197  3 

ZEROIN,  a  double  precision  FORTRAN  IV  subroutine 
to  solve  a  system  of  simultaneous,  nonlinear  equa- 
tions 

ZEROIN  computes  a  solution  vector  X  of  a  set  of  N 
simultaneous,  nonlinear  equations: 

Fj^(X)=0, 


METHOD: 


USAGE : 


F^(X)=0, 

X  is  obtained  by  an  iterative  method: 
j,(k+l)  ^  ^(k)^  ^(k) 

where  the  vector  D^'^'  =  J^'^^  *F(X^^' )  ,  with 
(k) 

(k)        (k) 
values  at  the  point  X    and  J    representing 


F(X    )  representing  the  vector  of  function 

)resi 

(k) 
an  approximation  to  the  Jacobian  at  X    .   The 

following  test  is  used  to  determine  convergence: 


||j^(k+l)  _  j^(k)  I  |2  ^  gpg^i  |j^ 


(k+1)  I  I  2 


where  EPS  is  a  user-supplied  tolerance.   All  real 
arithmetic  in  ZEROIN  is  double-precision. 

The  call  for  ZEROIM  is: 

CALL  ZEROIN (N,X,R, EPS, ICV,EVAL,IRITE) 

Input  values  are: 

N       the  order  of  the  system  of  equations, 

2  <  N  <  20. 

.  X       a  one-dimensional  array  initially  con- 
taining an  estimate  of  the  solution 
vector. 
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.  I 


AP-21-2 


Page  2 

R       an  estimate  for  | [solution  vector-x| | 

EPS     a  tolerance  to  be  used  in  the  test 
for  convergence,  usually  0.5E-12. 

IRITE    a  specification  for  intermediate 

printing.   If  IRITE  =  0,  ZEROIII  does 
no  printing.   If  IRITE  =  1,  ZEROIN 
prints,  at  each  iteration,  the  current 
value  of  R,  the  current  approximation 
to  the  solution,  the  number  of  times 
EVAL  (below)  has  been  called,  and  the 
current  function  values  F(X).   The 
printed  output  from  the  Tinal  itera- 
tion contains,  in  addition,  the  output 
value  of  ICV. 

EVAL    is  the  name  of  a  user-supplied 

SUBROUTItlE  subprogram  that  evaluates 
the  functions  F  at  each  iteration. 
The  call  for  EVAX.  is  CALL  EVAL{X,F) 
where  X  is  the  input  to  EVAL,  and  F 
contains  the  function  valu<bs  corres- 
ponding to  X.   X  and  F  roust  be  singly- 
dimensioned,  double-precision  vectors. 


Output  values  will  be: 


X 

R 


the  solution  vector. 


EPS 


ICV 


.(k+1) 


.(k) 


/Mx 


(k+1) 


where  X      is  the  solution  vector. 
When  ZEROIN  terminatoo  normally,  the  out- 
put value  of  R  will  be  less  than  the  in- 
itial value  of  EPS.   In  case  of  abnormal 
termination,  R  can  be  used  to  determine 
what  degree  of  convergence  was  obtained 
in  the  output  value  of  X. 

|If(X)||'  where  F(X)  is  the  vector  of 
function  values  corresponding  to  the 
solution  vector  X. 

one  of  the  following  indications: 

ICV  =0   if  convergence  was  obtained,  but 
I |F(X) I P  >  0.5E-10 


ICV 


ICV 


■  2 
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if  convergence  was  obtained,  and 
I l£(^) I  r  1  0.5E-10 

if  convergence  was  not  obtained 
within  M  iterations,  where 
M  =  MIN(10*N,50) . 


AP-21-2 


Appropriate  declarations  are; 


REAL* 8 
INTEGER* 4 
EXTERNAL 


XCN) ,  R,  EPS,  EVAL 
N,  ICV,  IRITE 
EVAL 


CODING  1. 

CONSIDERATIONS: 


AVAILABILITY: 


CORE 

REQUIREMENT: 

REFERENCES: 


AUTHOR: 


3. 


The  input  value  of  R  may  be  chosen  to  be 
I |x| 1/10.0,  where  X  is  the  input  vector  of  es- 
timates, if  the  solutions  are  all  about  the 
same  magnitude. 

The  output  valuiS  of  R  can  be  used  to  determine 
degree  of  convergence  in  the  case  ICV=2. 

ZEROIN  is  more  likely  to  converge  if  the  equa- 
tions are  scaled  so  all  solutions  are  of  the 
same  magnitude,  and  if  the  initial  estimates 
are  good.   Poor  initial  estimates  may  produce 
divide  check  interrupts  during  execution  of 
ZEROIN.   For  details  on  the  choice  of  initial 
estimates,  refer  to  the  first  reference  below. 


ZEROIN  is  available  in  load  module  form  in  the  li- 
brary SyS5.MATHLI3.SUDR  and  in  source  form  in 
SYS 5. J'lATHLIB. SOURCE.   For  further  information  about 
the  use  of  these  libraries,  refer  to  AP-60  or  con- 
sult the  Programming  Assistants  in  Room  39-219. 

On  the  370/165  ZEROIN  uses  9,730  bytes  of  core  stor- 
age (hexadecimal  2602) . 

Mancino,  O.G.  "Resolution  by  Iteration  of  Some  Non- 
linear Systems,"  Journal  of  the  Association  for 
Computing  Machinery,  Vol.  14,  No.  2  (April,  1967)  , 
pages  341-350. 

Robinson,  S.M.  "Interpolative  Solution  of  Systems 
of  Nonlinear  Equations ,  "  SIAiM  Journal  on  Numerical 
Analysis,  Vol.  3,  No.  4  (December,  1966),  page  650. 

ZEROIN  was  written  by  G.W.  Westley,  Computing  Tech- 
nology Center,  Union  Carbide  Corporation,  Nuclear 
Division,  Oak  Ridge,  Tennessee. 
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Computer  Flov/  Chart 
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Computer  Source  Listing 
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